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ABSTRACT 

We report the results from a systematic search for molecular (OH 119 fi m) outflows with Herschel/ PACS in a sample 
of 43 nearby (z < 0.3) galaxy mergers, mostly ultraluminous infrared galaxies (ULIRGs) and QSOs. We find that 
the character of the OH feature (strength of the absorption relative to the emission) correlates with that of the 9.7 //m 
silicate feature, a measure of obscuration in ULIRGs. Unambiguous evidence for molecular outflows, based on the 
detection of OH absorption profiles with median velocities more blueshifted than —50 km s _1 , is seen in 26 (70%) of 
the 37 OH-detected targets, suggesting a wide-angle (~145°) outflow geometry. Conversely, unambiguous evidence 
for molecular inflows, based on the detection of OH absorption profiles with median velocities more redshifted 
than +50 km s , is seen in only four objects, suggesting a planar or filamentary geometry for the inflowing gas. 
Terminal outflow velocities of 1000 km s _1 are measured in several objects, but median outflow velocities 
are typically 200 km s _1 . While the outflow velocities show no statistically significant dependence on the 
star formation rate, they are distinctly more blueshifted among systems with large active galactic nucleus (AGN) 
fractions and luminosities [log (Lagn/Tq) ^ 11.8 ± 0.3]. The quasars in these systems play a dominant role in 
driving the molecular outflows. However, the most AGN dominated systems, where OH is seen purely in emission, 
show relatively modest OH line widths, despite their large AGN luminosities, perhaps indicating that molecular 
outflows subside once the quasar has cleared a path through the obscuring material. 

Key words: galaxies: active - galaxies: evolution - ISM: jets and outflows - ISM: molecules - quasars: general 
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1. INTRODUCTION 

Gas-rich galaxy merging may trigger major starbursts, lead to 
the formation of elliptical galaxies, and account for the growth of 
supermassive black holes (BH; e.g., Sanders et al. 1988; Hopkins 
et al. 2009). This merger-driven evolutionary scenario starts with 
a completely obscured ultraluminous infrared galaxy (ULIRG). 
As the system evolves, the obscuring gas and dust is gradually 
dispersed, giving rise to dusty QSOs and finally to completely 
exposed QSOs. Powerful winds, driven by the central quasar 
or the surrounding starburst, have been invoked to stop the 
growth of both the BH and spheroidal component and explain 
the tight BH-spheroid mass relation (e.g., Fabian 1999; King 
2003; Murray et al. 2005). These winds are purported to quench 
star formation in the merger remnants (“negative mechanical 
feedback”), creating a population of red gas-poor ellipticals 
and explaining the bimodal color distribution observed in 
galaxy surveys (e.g., Kauffmann et al. 2003). There is growing 
observational support for these winds: e.g., most galaxies with 
high star formation rate (SFR) densities show signatures of 
outflows, both locally and at high redshifts (e.g., Heckman 2002; 
Veilleux et al. 2005; Chen et al. 2010; Weiner et al. 2009; Rubin 
et al. 2010; Steidel et al. 2010; Komei et al. 2012; Martin et al. 


2012; Newman et al. 2012, and references therein). The local 
outflows are often spatially resolved, allowing determination of 
the mass outflow rates (~0. 1—5 x SFR) and kinetic energies 
(~10 56 ~ 58 erg in ULIRGs; Rupke et al. 2002, 2005a, 2005b, 
2005c; Martin 2005, 2006). These winds are primarily driven 
by the starburst rather than the active galactic nucleus (AGN), 
except in late-stage mergers with quasar signatures where 
velocities >1000 km s _I (cf. 100-400 km s” 1 in other systems) 
are sometimes observed (e.g., Rupke et al. 2005c; Krug et al. 
2010; H. B. Krug et al. 2013, in preparation). 

Studies conducted in the past three years are shedding 
new light on the wind phenomenon in the local universe: (1) 
Early results from our Herschel 13 guaranteed time key program 
SHINING (PI: E. Sturm) have revealed far-infrared (FIR) OH 
features with P Cygni profiles indicative of massive molecular 
outflows in a number of ULIRGs (Fischer et al. 2010; Sturm 
et al. 201 1, hereafter F10 and SI 1, respectively; also Gonzalez- 
Alfonso et al. 2012, 2013). In a few of these objects, particularly 
those with dominant AGNs such as Mrk 231, the nearest quasar 


13 Herschel is an ESA space observatory with science instruments provided by 
European-led Principal Investigator consortia and with important participation 
from NASA. 
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known, the terminal outflow velocities exceed 1000 km s' 1 
and the outflow rates (up to ~1OOOM 0 yr” 1 ) are several times 
larger than the infrared (IR) based SFRs. (2) Independent, 
spatially resolved CO-emission observations of Mrk 23 1 with 
the IRAM/PdB millimeter- wave interferometer (Feruglio et al. 
2010; Cicone et al. 2012; it has also been mapped in HCN, 
FICO + , and HNC by Aalto et al. 2012) have confirmed this 
outflow and deduced mass outflow rates of ~700 M 0 yr -1 , far 
larger than the ongoing IR-based SFR (~16OM 0 yr -1 ) in the 
host galaxy. Further millimeter CO-observations have now been 
done by our group in a number of the objects discussed here, 
finding similar CO outflows (Cicone et al. 2013). Remarkably, 
the CO outflow in Mrk 231 partially overlaps spatially and 
kinematically with blueshifted optical Nai D Lk5890, 5896 
absorption features detected out to ~2-3 kpc from the nucleus 
(Rupke et al. 2005c). (3) Our high-resolution Gemini /Integral 
Field Unit (IFU) observations have revealed that the neutral- 
atom Na i D outflow in Mrk 23 1 is wide-angle, thus driven by a 
QSO wind rather than a jet (Rupke & Veilleux 2011, hereafter 
RV11; for short, neutral-atomic Nai D outflows are simply 
called neutral outflows in the rest of this paper, unless noted 
otherwise). Similar wide-angle high-velocity outflows have now 
been mapped at high resolution in a number of other local 
ULIRGs (Rupke & Veilleux 2013a, hereafter RV13, 2013b). 

These powerful outflows may be the long-sought “smoking 
gun” of quasar mechanical feedback that clears out the molec- 
ular disk formed from dissipative collapse during the merger. 
The tentative trend of increasing OFI terminal outflow velocity 
with increasing AGN luminosity identified by SI 1 adds support 
to this idea, but it is based on only six ULIRGs (and NGC 253, 
a regular starburst galaxy). To properly test this idea, one needs 
to examine a statistically representative sample of ULIRGs and 
study the properties of their molecular outflows uniformly and 
systematically. In this paper, we report the first results from 
such a study, focusing on the directly measured outflow ve- 
locities (in a later paper, we will report on our analysis of the 
model-dependent dynamical quantities (masses, momenta, and 
energies) in the subset of objects for which we have multiple 
OH transitions). The sample used for this analysis is described 
in Section 2. The observations and methods used to reduce 
and analyze the data are discussed in Section 3. The results 
from this analysis are presented in Section 4 and their impli- 
cations are discussed in Section 5. The main conclusions are 
summarized in Section 6. Throughout this paper, we adopt Hq — 
70 km s -1 Mpc" 1 , £l M = 0.3, and £2 A = 0.7. We also adopt the 
standard convention that approaching material has a negative 
velocity with respect to the systemic velocity of the host galaxy 
due to Doppler shift. 

2. SAMPLE 
2.1. Sample Selection 

As mentioned in Section 1 , the first results from the SHINING 
GTO survey of ULIRGs (S 1 1 ) revealed a tentative trend between 
OH outflow velocities and AGN luminosities, but this was 
largely driven by the velocities measured in Mrk 231 and 
F08572+3915, the only two warm (large / 25//60 flux ratio), 
quasar-dominated, late-stage mergers in the Sll sample. It 
quickly became clear that the rest of the SHINING GTO data 
would not be able to address this question adequately because 
the GTO sample largely probes the early stages of the merger 
after the first peri-passage: 20 of the 23 (87%) GTO ULIRGs are 
cold (small //s //$ 0 ratio) pre-mergers or recent mergers, where 


the bolometric and mechanical luminosities are still dominated 
by the nuclear starburst (Veilleux et al. 2009b, hereafter V09). A 
cycle 1 open-time program (OTl_sveilleu_l, PI: S. Veilleux) has 
allowed us to expand the original GTO sample to include a set of 
warm quasar-dominated ULIRGs. These objects were selected 
using four criteria. (1) The targets had to be part of the QUEST 
sample of local (z. < 0.3) luminous mergers: either ULIRGs 
or QSOs (Veilleux 2012, and references therein). (2) Their 
bolometric luminosity had to be dominated by the quasar based 
on the Spitzer data or, equivalently, //M.S'-band 25 to 60 pm flux 
ratio / 25//60 Pi 0-15 (this criterion also automatically selects 
late-stage, fully coalesced mergers; see V09). (3) High signal- 
to-noise ratio (S/N ~ 30-40) Nai D absorption spectra had to 
be planned or available for each object. This criterion did not 
bias the sample in any way since these objects were selected 
randomly based on visibility at the time of the ground-based 
observations. (4) They had to have /100 > 1 Jy so that high S/N 
in the continuum could be reached in a reasonable amount of 
time with Herschel- PACS. These criteria resulted in a sample of 
18 objects, 3 of which were already part of the SHINING GTO 
survey of ULIRGs. This sample size was a good match to the 
existing 20 cold pre-merger ULIRGs in SHINING. 

Since all ULIRGs in the SHINING survey and its OT1 
extension were selected to have / AAV band 100 pm fluxes 
/100 > 1 ly to more easily probe OH 119/xm, this sample 
did not contain any “classic” IR-faint QSOs in the critical 
late merger phases when the quasar has finally gotten rid of 
its natal “cocoon” and the effects of feedback are predicted 
to subside (Narayanan et al. 2008; Hopkins et al. 2009). 
Therefore we put in a successful open-time request in Cycle 2 
(OT2_sveilleu_4, PI: S. Veilleux) to further expand the sample by 
including objects meeting the following three additional criteria. 

(1) IR-to-bolometric luminosity ratios smaller than 0.6 i.e., 
IR-fainter than all ULIRGs in the SHINING and OT1 programs 
(IR here and throughout this paper refers to 8-1000 pm). 

(2) High-S/N (~30— 40) Nai D and FUV absorption spectra 
had to be planned or available for each object. The new FUV 
criterion favored the “classic” UV-bright QSOs. (3) The /100 
threshold was lowered by a factor of two down to >0.5 Jy to 
accommodate these IR-fainter systems. This factor of two was 
chosen to provide enough objects in the OT2 sample without 
requiring unrealistically long exposure times per target. These 
criteria resulted in a sample of five objects. While this sample 
size is small, recall that the main objective of this OT2 extension 
is to anchor the merger sequence at the latest (IR-faintest) stages 
that could be probed with Herschel. This is a merger phase that 
has not been explored with Herschel until now so each of these 
objects provides important new information. 

In summary, the sample consists of 43 objects, 23 from the 
GTO sample, 15 from the OT1 program, and 5 from the OT2 
program. All of these objects show signs of ongoing or recent 
interactions (NGC 4418 = F12243— 0036, a GTO target, is not a 
major merger but it is IR luminous so we include it in the present 
sample). We include here the ULIRGs from SI 1 again, applying 
the same fitting method as for the rest of the sample, in order to 
be complete and consistent. The main objective of the present 
study is to search for trends with AGN and host properties as 
well as age across the merger sequence (after the first encounter, 
post-LIRG phase). In Section 3, we will subdivide these 43 
mergers into various subcategories (e.g., morphology, starburst 
power, AGN luminosity) to try to find the driving parameter of 
the wind properties. Our Spitzer results (V09) have shown that 
BH accretion, in addition to depending on the merger phase. 
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Figure 1 . Histograms showing the distributions of the sample properties: (a) redshifts, (b) starburst luminosities, (c) AGN fractions, and (d) AGN luminosities. 


also has a strong chaotic/random component, so a sample size 
substantially smaller than 43 objects would make it hard to draw 
statistically significant conclusions. 

2.2. Properties of the Sample Galaxies 

The properties of our sample galaxies are listed in Table 1. 
The notes to Table 1 briefly explain the meaning of each of these 
quantities. The luminosities listed in Table 1 rely on a number of 
assumptions that deserve further explanation. Following V09, 
the bolometric luminosities for ULIRGs were estimated to be 
CbOL = 1.15 Ljr , where L IR is the IR luminosity over 8-1000 p m 
(Sanders & Mirabel 1996), and Lbol = 7L(5100 A) + Lir 
for the Palomar-Green (PG) QSOs (Netzer et al. 2007). The 
starburst and AGN luminosities were next calculated from 

Ebol = 7-agn + Esb (1) 

= “agnEbol + Esb, (2) 

where q!agn is the fractional contribution of the AGN to the bolo- 
metric luminosity, hereafter called the “AGN fraction” for short. 
In V09, we measured o!agn in the QUEST ULIRGs and QSOs 
using six independent methods that span a range in wavelength 
and give consistent results within ~±10%-15% on average. 
These six methods relied on (1) the [Oiv]/[Nen] ratio, (2) the 
[Nev]/[Neii] ratio, (3) the polycyclic aromatic hydrocarbon 
(PAH) 7.7 pm equivalent width, (4) the modified version of 
the Laurent et al. (2000) diagram which plots the 14-16 pm to 
5. 3-5 . 8 pm continuum ratio versus the PAH 6.2 pm equivalent 
width, (5) the PAH- free, silicate-free 5-25 pm to FIR continuum 
ratio, and ( 6 ) the 15 to 30 pm continuum ratio, / 15 // 30 . 

Unfortunately, several objects in the SHINING GTO sample 
are not part of the QUEST sample and therefore do not have 
such accurately measured AGN contribution to the bolometric 
luminosity. In the following discussion, we rely solely on 
the restframe / 15//30 continuum ratio, which is available for 
all the objects in the Herschel sample, to quantify the AGN 
contributions to the bolometric luminosities of our systems. 


This ratio was found by V09 to be more tightly correlated 
with the PAH-free, silicate-free mid-IR (MIR) /FIR ratio and the 
AGN contribution to the bolometric luminosity than any other 
Spitzer - derived continuum ratio. 

Following V09, we use / 30//15 as a surrogate of the 
PAH-free, silicate-free MIR/FIR ratio and adopt log(/ 3 o//i 5 ) = 
0.20 and 1.35 as the zero points for the “pure” AGN and the 
“pure” starburst ULIRG, respectively (note that V09 make a dis- 
tinction between “pure” starburst ULIRGs and “normal” star- 
forming galaxies). The AGN contribution is calculated from a 
linear interpolation between these two extremes. The zero point 
for the pure AGN corresponds to the average / 30//15 ratio of 
FIR-undetected PG QSOs, while the zero point for the pure 
starbursts is calculated from the 10 ULIRGs with the largest 
/ 30//15 ratios. We used the correction factors listed in Table 10 
of V09 to transform 15 pm luminosities into bolometric lu- 
minosities: log[E v (15 /xm)/Z.BOL] = —14.33 for a pure AGN 
and —14.56 for a pure starburst ULIRG. The uncertainty on 
“agn is estimated to be ± 20 % on average, but is likely higher 
for strongly buried sources with significant FIR extinctions 
(e.g., Arp 220 = F15327+2340). We return to the issue of FIR 
extinction in Sections 4.6 and 5. 

Figure 1 shows the distributions of redshifts, starburst 
luminosities, AGN fractions, and AGN luminosities for all 43 
objects in our sample. This figure emphasizes the broad range of 
properties of our sample galaxies. The starburst and AGN lumi- 
nosities are independent of each other i.e., we find no correlation 
or anti-correlation between the starburst and AGN luminosities 
in our sample of galaxies, but perhaps not surprisingly objects 
with the largest starburst (AGN) luminosities (>10 12 L 0 ) tend 
to have lower (higher) AGN fractions. 

3. OBSERVATIONS, DATA REDUCTION, 

AND SPECTRAL ANALYSIS 

3.1. Observations 

All of the data in this paper were obtained with the PACS 
FIR spectrometer (Poglitsch et al. 2010) on board Herschel 
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Name 

Other Name 

z 

111 

“AGN 

log Lbol 

log isB 

log Eagn 

<7* 

M h 

IC 

Type 

Program (/ e xp) 





(%) 

(Lq) 

(Lq) 

(Lq) 

(km s -1 ) 

(mag) 



(hr) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ID 

(12) 

(13) 

F00509+1225 

IZw 1 

0.0610 

0.434 

90.1 

12.01 

11.00 

11.96 

188 


s 

AGN 1 

OT1 (1.3) 

F01572+0009 

Mrk 1014 

0.1631 

0.161 

64.6 

12.68 

12.23 

12.49 

200 

-25.39 

s 

AGN 1 

OT1 (1.0) 

F05024— 1941 


0.1920 

0.050 

7.3 

12.43 

12.40 

11.30 


-24.49 

B 

AGN 2 

OT1 (4.6) 

F05 189-2524 


0.0426 

0.198 

71.7 

12.22 

11.67 

12.07 

137 

-23.96 

s 

AGN 2 

GTO (0.5) 

07251-0248 


0.0876 

0.075 

30.0 

12.45 

12.29 

11.92 



s 


GTO (1.1) 

F075 99+65 08 


0.1483 

0.377 

87.6 

12.59 

11.68 

12.53 


-24.55 

s 

AGN 1 

OT1 (3.7) 

F08572+3915 


0.0584 

0.191 

70.4 

12.20 

11.67 

12.05 


-23.58 

B 

L 

GTO (1.0) 

09022-3615 


0.0596 

0.125 

54.9 

12.35 

12.00 

12.09 



s 


GTO (0.5) 

F09320+6134 

UGC5101 

0.0394 

0.130 

56.4 

12.05 

11.69 

11.80 



s 

L 

GTO (0.3) 

F10565+2448 


0.0431 

0.105 

47.1 

12.10 

11.82 

11.77 



B 

H 

GTO (0.5) 

FI 11 19+3257 


0.1890 

0.266 

79.9 

12.71 

12.02 

12.61 


-23.44 

s 

AGN 1 

OT1 (4.6) 

Z1 1598-01 14 


0.1507 

0.210 

73.5 

12.60 

12.03 

12.47 


-25.06 

s 

AGN 1 

OT1 (0.8) 

FI 2072-0444 


0.1286 3 

0.220 

74.8 

12.46 

11.86 

12.33 


-24.22 

B 

AGN 2 

OT1 (0.8) 

F121 12+0305 


0.0733 

0.060 

17.8 

12.38 

12.30 

11.63 


-24.36 

B 

L 

GTO (0.5) 

F12243— 0036 

NGC4418 

0.0071 

0.128 

55.7 

11.06 

10.71 

10.81 



S 

AGN 2 

GTO (0.2) 

F12265+0219 

3C 273 

0.1583 

0.840 

98.5 

12.86 

11.04 

12.86 


-28.06 

s 

AGN 1 

OT1 (1.0) 

FI 2540+5708 

Mrk 231 

0.0422 

0.272 

80.5 

12.60 

11.89 

12.51 

120 

-24.52 

s 

AGN 1 

GTO (0.3) 

13120-5453 


0.0308 

0.080 

33.4 

12.30 

12.12 

11.83 



s 

AGN 2 

GTO (0.3) 

F13305— 1739 


0.1484 

0.377 

87.5 

12.32 

11.42 

12.26 


-25.47 

s 

AGN 2 

OT1 (4.9) 

F13428+5608 

Mrk 273 

0.0378 

0.081 

34.2 

12.21 

12.03 

11.74 

285 

-24.32 

s 

AGN 2 

GTO (0.3) 

F13451+1232 


0.1217 

0.273 

80.6 

12.38 

11.67 

12.29 

146. 167 b 

-25.52 

B 

AGN 2 

OT1 (1.0) 

F14348— 1447 


0.0830 

0.060 

17.4 

12.40 

12.31 

11.64 

160 

-24.99 

B 

L 

GTO (1.1) 

F14378— 3651 


0.0676 

0.064 

21.1 

12.17 

12.07 

11.50 

153 


s 

AGN 2 

GTO (1.1) 

F14394+5332 


0.1045 

0.152 

62.5 

12.17 

11.74 

11.96 


-24.88. — 23.56 c 

B b 

AGN 2 

OT1 (1.0) 

FI 5206+3342 


0.1244 

0.158 

63.9 

12.30 

11.86 

12.10 


-24.47 

S 

H 

OT1 (1.5) 

FI 5250+3608 


0.0552 

0.095 

42.3 

12.10 

11.86 

11.73 

150 


S 

L 

GTO (0.5) 

F15327+2340 

Arp 220 

0.0181 

0.049 

5.8 

12.22 

12.19 

10.98 

164 


B 

L 

GTO (0.9) 

F15462— 0450 


0.0998 

0.145 

60.6 

12.27 

11.86 

12.05 

169 

-23.4 

S 

AGN 1 

OT1 (1.0) 

FI 6504+0228 

NGC 6240 

0.0245 

0.130 

56.2 

11.90 

11.54 

11.65 

229 


B 

L 

GTO (0.3) 

F17207— 0014 


0.0428 

0.038 

<5% 

12.45 

12.43 

<11.15 

229 


S 

L 

GTO (0.3) 

F19297— 0406 


0.0857 

0.066 

23.4 

12.44 

12.33 

11.81 



S 

H 

GTO (0.3) 

19542+1110 


0.0624 d 

0.069 

25.5 

12.12 

11.99 

11.52 



s 

L 

GTO (1.1) 

F2055 1-4250 


0.0430 

0.132 

56.9 

12.11 

11.74 

11.87 

140 


s 

L 

GTO (0.5) 

F2249 1 — 1808 


0.0778 

0.057 

14.4 

12.90 

12.83 

12.05 


-24.21 

B 

H 

GTO (0.8) 

F23128— 5919 

ESO 148-IG 002 

0.0446 

0.154 

63.0 

12.09 

11.66 

11.89 



B 

H 

GTO (0.3) 

F23233+2817 


0.1140 

0.204 

72.5 

12.11 

11.55 

11.97 


-24.54 

S 

AGN 2 

OT1 (1.0) 

F23365+3604 


0.0645 

0.100 

44.6 

12.22 

11.97 

11.87 

145 


s 

L 

GTO (0.5) 

F23389+0300 


0.1450 

0.066 

22.7 

12.19 

12.07 

11.54 


-24.15 

B 

AGN 2 

OT1 (4.4) 

PG 1126-041 

Mrk 1298 

0.0620 

0.329 

88.7 

11.52 

10.57 

11.47 

194 

-24.15 

s 

AGN 1 

OT2 (3.1) 

PG 1351+640 


0.0882 

0.389 

92.2 

12.04 

10.93 

12.00 



s 

AGN 1 

OT2 (10.6) 

PG 1440+356 

Mrk 0478 

0.0780 

0.361 

90.7 

11.80 

10.77 

11.76 


-24.65 

s 

AGN 1 

OT2 (5.4) 

PG 1613+658 

Mrk 0876 

0.1290 

0.398 

92.6 

12.29 

11.16 

12.26 


-25.68 

s 

AGN 1 

OT2 (4.1) 

PG 2130+099 

UGC 11763 

0.0630 

0.563 

98.5 

11.77 

9.95 

11.76 

156 

-24.14 

s 

AGN 1 

OT2 (14.0) 


Notes. Column 1: galaxy name. Coordinate-based names beginning with F are sources in the IRAS Faint Source Catalog. Column 2: another name. Column 3: redshift. 
Column 4: /15//30 values from V09 and Stierwalt et al. (2013). When no 15 to 30 pm data is available the value is estimated by using the /25//6O ratio and the good 
correlation between the 25 to 60 versus 30 to 15 pm flux ratios (see V09). Column 5: c^agn^ fractional contribution of the AGN to the bolometric luminosity based on 
the /i 5 // 3 o- See Section 2.2. Column 6: bolometric luminosity. See Section 2.2 for more detail. Column 7: starburst bolometric luminosity. See Section 2.2 for more 
detail. Column 7: AGN luminosity. See Section 2.2 for more detail. Column 8: o *, stellar velocity dispersion from Dasyra et al. (2006a, 2006b, 2007). Column 10: 
Mh, absolute H - band magnitudes from Veilleux et al. (2006, 2009b). For the sake of completeness we also included Mh absolute magnitudes derived from M p 
values assuming H-K ' = 0.35 mag (Veilleux et al. 2002). Column 1 1 : interaction class, “S” for singles and “B” for binaries, based on classifications of Veilleux et al. 
(2002, 2006, 2009b) and, when not available, our own visual classification from //5T/ACS F814W images (Iwasawa et al. 201 1). Column 12: optical spectral types 
of our sample, Hu (for H 11 galaxies), L (for LINER-like), AGN 2 and AGN 1 (for type 2 AGN and type 1 AGN, respectively). For the spectral type we adopted the 
classification from Veilleux et al. (1999) and V09, when not available, we used the values from NED/SIMBAD. Column 13: program, GTO SHINING program (PI: 
E. Sturm), QUEST OT1 and OT2 programs (PI: S. Veilleux). In parenthesis is the total PACS exposure time. 

a Redshift from Spoon & Holt (2009). This value is based on [Ne 11 ] 12.8 pm and is slightly larger than that derived from optical data (e.g., Strauss et al. 1992). 
b The values presented correspond to the eastern and western nucleus of FI 345 1+1232, respectively. 
c F14394+5332 is a triple system with an eastern pair and a western component. 

d Average redshift derived from our fit to the [Neil], [Nem] and [Oiv] emission lines in the Spitzer / IRS spectra. 
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(Pilbratt et al. 2010). Several systems in the OT1 and OT2 
sub-samples are relatively faint in the FIR so an exhaustive 
multi-line PACS survey was not possible. Contrary to the GTO 
data, which generally cover the ground-state OH 119 /tm 2 Fl 3/2 
J — 5/2— 3/2 rotational A-doublet transitions, the high-lying 
(Slower = 290 K) 65 /x m 2 IT 3/2 J —1 /2— 9/2 and (£i OW er = 120 K) 
84 fim 2 n 3/2 J =1 /2—S/2 rotational A-doublet transitions, and 
the ground-state cross-ladder 79 /xm 2 ni/2- 2 Il3/2 J = 1/2-3/2 
rotational A-doublet transitions, the OT1 and OT2 data focus on 
a single OH feature, the OH 1 19 /xm doublet (and ls OH 120 /im 
counterparts). The value of this feature as a wind diagnostic 
was beautifully demonstrated in F10 and Sll and our early 
GTO data. It is the strongest OH transition in Mrk 231 and 
most of the GTO targets. It is positioned in wavelength near 
the peak spectroscopic sensitivity of PACS. Finally, it is seen 
in absorption more often than the 79 pm feature and thus is 
more likely to provide an unambiguous signature of gas motion 
(rotation, inflow, or outflow). 

PACS was used in range scan spectroscopy mode in high 
sampling centered on the redshifted OH 119/xm + 18 OH 
120 /xm complex with a velocity range of ~8000 km s' 1 
(restframe 118-121 pm ) to provide enough coverage on both 
sides of the OH complex for reliable continuum placement. The 
resulting PACS spectral resolution is ~270 km s -1 . The total 
amount of time spent to carry out these observations (including 
overheads) was 13.0, 32.5, and 37.2 hr during GT, OT1, and 
OT2, respectively. The on-target exposure times were adjusted 
according to the continuum level under OH, estimated from an 
extrapolation of the Spitzer MIR + IRAS FIR energy distribution. 
A small chopper throw of 1 .'5 was used in most cases. The 
observing time for each target, including all overheads, is listed 
in the last column of Table 1 . 

3.2. Data Reduction 

All of the PACS data were reduced in the same way as in 
Sll. As described there, the data reduction was done using 
the standard PACS reduction and calibration pipeline (ipipe) 
included in HIPE 6.0. The spectra were normalized to the tele- 
scope flux (which dominates the total signal) and recalibrated 
it with a reference telescope spectrum obtained from dedicated 
Neptune observations during the Herschel performance verifi- 
cation phase. All of our objects are point sources for PACS. In 
the following, we use the spectrum of the central 9" x 9" spatial 
pixel (spaxel) only, applying the point-source correction factors 
(PSF losses) as given in the PACS documentation. We have ver- 
ified this approach by comparing the resulting continuum flux 
density level to the continuum level of all 25 spaxels combined 
(which is free of PSF losses and pointing uncertainties). In all 
cases the agreement is excellent, however the central spaxel 
alone provides better S/N. 

The reduced spectra were next smoothed using a Gaussian 
kernel of width 0.05 pm (i.e., about half a resolution element) 
to reduce the noise in the data before the spectral analysis. 
A spline was fit to the continuum and subtracted from the 
spectra; these continuum-subtracted spectra were used for the 
subsequent spectral fitting. 

3.3. Spectral Analysis 

The profiles of the OH 119.233, 119.441 pm doublet were 
modeled using four Gaussian components. The fits were carried 
out using PySpecKit, an open-source, extensible spectroscopic 
analysis toolkit for astronomy (Ginsburg & Mirocha 2011). 


This tool uses the Levenberg-Marquardt technique to solve 
the least-squares problem in order to find the best fit for the 
observations. 14 Each line of the OH doublet was fitted with two 
Gaussian components characterized by their amplitude (either 
negative or positive), peak position, and standard deviation (or, 
equivalently, FWHM). The separation between the two lines of 
the doublet was set to 0.208 p m in the restframe (~520 km s -1 ) 
and the amplitude and standard deviation were fixed to be 
the same for each component in the doublet. Overlap effects 
between the various components were ignored for simplicity 
(e.g., the emission component in P Cygni profiles is not affected 
by the absorption component in the foreground. This is discussed 
further in the next paragraph). The 18 OH 120/xm doublet is 
detected in a number of objects but is too far to the red to 
significantly affect the profiles of the OH 119/xm feature. 
Likewise, contamination from 17 OH, H20 + , CH, CH + , and CO 
(22-21) is considered unlikely since these lines are either too 
far to the red or blue of the OH 1 19 pm feature or undetected 
even in the high-S/N spectrum of Mrk 231 (see Figure 2(a) of 
Gonzalez-Alfonso et al. 2013). 

We find that three distinct scenarios apply to our data: (1) pure 
OH absorption, (2) pure OH emission, and (3) P Cygni profiles 
(our data reveal no new evidence for inverted P Cygni OH pro- 
files; this interesting result is discussed further in Section 4.3). In 
scenario 1, there is no evidence for any OH emission and each 
line of the OH doublet is fitted with two absorption compo- 
nents, tracing both the stronger low-velocity component of the 
outflow and the fainter high-velocity component (e.g., Mrk 23 1 ; 
Gonzalez-Alfonso et al. 2013). Scenario 2 is treated similarly. 
In this scenario, there is no evidence for any OH absorption 
feature and two Gaussian components are used to model each 
of the line in the OH doublet. Scenario 3 is more subtle. We first 
tried to fit the P Cygni profiles with three or four components per 
line (1-2 blueshifted absorption components and 1-2 redshifted 
emission components), but we quickly realized that these fits 
were underconstrained. We had to settle on fitting these profiles 
with a single blueshifted absorption component and a single red- 
shifted emission component i.e., no attempt is made to capture 
the fainter high-velocity component of the outflow. This has the 
potential to underestimate the outflow velocities derived from 
this method (making the velocities derived from the absorption 
feature more positive than they should be). There is also the risk 
of a degeneracy between the emission and absorption compo- 
nents since in principle they each could be unrealistically strong 
as long as they cancel each other to reproduce the spectrum. In 
the radiative transfer models of Gonzalez-Alfonso et al. (2012, 
2013), overlap effects between the background emission com- 
ponent and the foreground absorption component prevent this 
degeneracy, virtually zeroing out the emission component(s) 
over the range of velocities of the absorption component(s). To 
approximate these overlap effects, the fit in some cases is fur- 
ther constrained by requiring the absorption feature to not be 
deeper than the observed absorption. This assumption was made 
for six objects (F09022-3615, F10565+2448, F12072-0444, 
F121 12+0305, F14348— 1447, and F23389+0300), but in the 
end there are only two cases where the emission component 
affects the absorption from the blue (119.233 um) line of the 
doublet: F121 12+0305 and F23389+0300. However, even in 
these objects, the emission contribution is weak at the posi- 
tion of this line and does not affect the measured velocities by 


14 By default PySpecKit implements the Levenberg-Marquardt algorithm via 
MPFIT (Markwardt 2009). 
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more than the typical measurement uncertainties (±50 km s -1 ; 
see below), based on comparisons with unconstrained fits. In 
Section 4.3, we compare the velocities measured in scenario 2 
with those measured in scenario 1 and find no significant shift, 
suggesting that the assumptions made in scenario 2 have no 
systematic effect on the velocity measurements. Nevertheless, 
we use caution and distinguish the absorption velocities derived 
under scenario 2 from those measured in scenarios 1 and 3 in 
many of the figures presented in Section 4. 

These fits were first used to quantify the strength and nature 
(absorption versus emission) of this feature: (1) the total flux and 
equivalent width of the OH 1 19.441 /xm line, adding up all of the 
absorption and emission components, (2) the flux and equivalent 
width of the absorption component(s) used to fit this line, and 
(3) the flux and equivalent width of the emission component(s) 
used to fit this line. 

We next characterized the OH profile by measuring a few 
characteristic velocities from our fits: (1) i>5o(abs), the median 
velocity of the fitted absorption profile, i.e., 50% of the ab- 
sorption takes place at velocities above (more positive than) 
U5o(abs), (2) i>84(abs), the velocity above which 84% of the ab- 
sorption takes place, (3) i>5o(emi), the median velocity of the 
fitted emission profile, i.e., 50% of the emission takes place at 
velocities below (less positive than) uso(emi), and (4) U84(emi), 
the velocity below which 84% of the emission takes place. 
We did not extend our fit-based analysis beyond 84% of the 
absorption/emission profiles because the wings of the OH pro- 
files may not be well captured by our simple two-component fits 
(especially in the case of P Cygni profiles where a single Gaus- 
sian was used to model the absorption (or emission) feature of 
each line in the doublet). For comparison with SI 1, we also es- 
timated the terminal outflow velocity, i> max , from the maximum 
extent of the blueshifted wing of the OH 1 19.233 /xm absorption 
profile (it was not measured in the case of pure OH emission 
profiles). This quantity was measured independently by three 
members of our team to try to reduce possible biases, but it is 
admittedly more subjective and uncertain than the other velocity 
measurements, as noted below. 

There are a number of potential sources of uncertainties when 
measuring the equivalent widths, fluxes, and velocities from our 
data. The wavelength uncertainty depends on mispointing. Any 
mispointing may also have an effect on the line profiles. We 
checked for this and found the effect not to apply in our cases. 
We find instead that the uncertainties on these measurements are 
dominated by the placement of the continuum. To estimate this 
source of uncertainty, we re-measure these quantities assuming 
different continuum shapes (spline versus polynomial) and 
positions. We find that the typical uncertainties on the equivalent 
widths and fluxes are ±20% and ±50 km s _1 on the velocities, 
except for u max which is much more uncertain (±200 km s” 1 ). 
For this reason, we give considerably more weight to results 
based on vso and Vs 4 in the following discussion. 

4. RESULTS 

Figure 2 shows the fits to the OH 119/xm profiles. The 
OH 119 /xm parameters derived from these fits are tabulated in 
Table 2 along with their uncertainties. The average and median 
values of these parameters are listed at the end of this table. The 
meaning of each parameter is discussed in Section 3.3 and the 
notes to Table 2. Note that the fluxes and equivalent widths in that 
table need to be multiplied by a factor of two when considering 
both lines of the doublet. In this section we compare these results 
with the other galaxy properties tabulated in Table 1 . 


4.1. Detection Rate of the OH 119 /xm Feature 

The 119/xm OH doublet was detected in 37 of the 43 objects 
in our sample (86%). The OH detection rate shows no obvious 
dependence on the properties of the objects (Figure 3). A closer 
inspection of the data indicates that the six objects without OH 
detection have systematically lower S/N ratios (<20-30) in 
the FIR continuum, and therefore less stringent limits on the 
OH 119/xm equivalent widths, than the others. Two of them 
are IR-faint QSOs from the OT2 sample (PG 1126—041 and 
PG 2130+099). The other four are part of our OT1 sample of 
warm ULIRGs (F13305-1739, F15206+3342, F23128-5919, 
and F12265+0219 = 3C 273). These six objects are excluded 
from the rest of the analysis. 

4.2. Nature of the OH 119 /im and 
Mid-infrared Silicate Features 

Of the 37 OH 119/xm detections, 17 are seen purely in 
absorption, 15 15 show absorption + emission composite profiles, 
and 5 are purely in emission. Figure 4 shows no obvious 
trend between the total equivalent width of the OH feature 
and AGN fraction, starburst luminosity or AGN luminosity, 
except for the fact that all four objects with AGN fractions 
“agn > 90% show pure OH emission. This result is consistent 
with the detection of OH 119, 79, and 163 /xm emission in 
the prototypical Seyfert 2 galaxy NGC 1068 (Spinoglio et al. 
2005; Hailey-Dunsheath et al. 2012). Preliminary results on the 
Swift Burst Alert Telescope AGN sample (OT2 sveilleu 6, PI: 
S. Veilleux), where OH 1 19 /xm is often seen in emission but 
log (Lagn/±o) < 12 (M. Melendez et al. 2013, in preparation), 
reinforces the idea that the AGN fraction is more important 
in setting the character (i.e., strength of emission relative to 
absorption) of the OH feature than the AGN luminosity. These 
results suggest that dominant AGN in ULIRGs and QSOs 
provide a favorable environment to produce the OH molecule 
and excite it (via radiative pumping or collisional excitation) to 
the ground-state 2 n3/2 J = 5/2 level. Interestingly, a similar 
trend with AGN fraction was seen by Teng et al. (2013), 
when examining the Hi 21 cm feature in a smaller set of 
QUEST ULIRGs and QSOs with the Green Bank Telescope 
(GBT). 

Figure 5 compares the strength of the OH 119 /xm absorption 
feature with that of the 9.7 /xm silicate absorption feature, a 
measure of obscuration in ULIRGs. The strengths of the 9.7 /xm 
feature in our sample galaxies are taken from two separate papers 
that use two different methods to derive this quantity. In V09, 
T9.7 fim is the effective 9.7 /xm silicate optical depth measured 
relative to the sum of the fitted blackbody components (a larger 
value of t9.7 Mm implies a deeper silicate absorption feature). 
In Stierwalt et al. (2013), Sg -j ^ is the logarithm of the ratio 
of the measured flux at the central wavelength of the silicate 
absorption feature to the local continuum flux (a more negative 
value of Sgj^m implies a deeper silicate absorption feature). 
In both cases, we note that the more obscured systems (T9.7 
> 4 in V09 or S9.7 ;ml < —2 in Stierwalt et al. 2013) more 
often show OH purely in absorption than the less obscured 
systems. Conversely, systems with strong OH emission either 
show very weak silicate absorption or silicate in emission. This 
last category includes three of the four QSOs seen in emission 
in both OH and silicate (F00509+1225 = IZw 1,PG 1440+356, 


15 Some of these 17 objects (e.g., F05024— 1941, F15462— 0450, and 
F2055 1—4250) may have faint redshifted emission but our simple fitting 
procedure was not successful in capturing this emission. 
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Figure 2. Spectral fits to the OH 119 (tm profiles in the 43 objects of our sample. In each panel, the solid black line represents the data and the solid purple line is 
the best multi-component Gaussian fit to these data. The blue dash line represents the absorption component(s) used in this fit while the red dash line represents the 
emission component(s). The origin of the velocity scale corresponds to OH 1 19.233 (tm at the systemic velocity. The two vertical dashed and dotted lines mark the 
positions of the 16 OH and I8 OH doublets, respectively. The vertical dot-dash line marks the position of CH + 1 19.848 //m. 

(A color version of this figure is available in the online journal.) 


and PG 1 6 13+658; PG 1351 +640 is the only QSO that is not part 
of the Spitzer sample of Schweitzer et al. 2008; these objects are 
not plotted in Figure 5, but would reinforce the trend between 
FIR OF! and MIR silicate). 

Figures 4 and 5 paint a consistent picture where both the 
OH 1 19 pm and 9.7 pm silicate features change character from 
absorption to emission as the merger progresses and the AGN 
becomes dominant. 


4.3. Distributions of Velocities and Wind Detection Rates 

Figure 6 shows the distributions of velocities derived from 
both the OH absorption and emission line features (uso(abs), 
i> 84 (abs), i) 5 o(emi), i> 84 (emi), and u max (abs), as defined in 
Section 3.3 and listed in Table 2). The velocity distributions 
of the absorption features show a distinct excess at negative 
values while the opposite is true of the emission features. As 
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mentioned in Section 3.3, all 15 objects with clear composite 
absorption + emission features show a distinct P Cygni profile 
i.e., a blueshifted absorption feature accompanied with a red- 
shifted emission feature, an unambiguous signature of outflow. 
As listed at the bottom of Table 2, the average (median) loofahs), 
U 84 (abs), and i> max (abs) are —194 (—204), —444 (—492), and 
—927 (—925) km s -1 , respectively. Given possible projection 
effects, these measurements represent lower limits on the true 
velocities of the OH outflows in these objects. 

Figures 6(c) and (d) show i> 5 o(abs) and i> 84 (abs) in objects 
with pure absorption features and P Cygni profiles. The aver- 
age (median) i> 5 o(abs) and i> 84 (abs) are —135 (—183) and —398 


(—477) km s~ 1 for the objects with pure OH absorption profiles 
versus —262 (—243) and —495 (—495) km s -1 for the objects 
with P Cygni profiles. Given the measurement uncertainties 
(±50 km s -1 or larger in some cases; Table 2), we do not 
consider these differences significant. A Kolmogorov-Smirnov 
(K-S) between these distributions supports this conclusion. As 
a case in point, we generally find excellent agreement between 
our velocities and those of Sll and Gonzalez- Alfonso et al. 
(2012; the largest discrepancy is seen for the P Cygni profile of 
F14387— 3651, where Us 4 in SI 1 is 240 km s -1 more blueshifted 
than the value reported here). Nevertheless, to remain aware 
of possible biases, the velocities derived in objects with 
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OH P Cygni profiles will be labeled differently from those de- 
rived in objects with OH purely in absorption or emission in all 
of the following figures. 

We follow Rupke et al. (2005b) and conservatively define a 
wind as having an OH absorption feature with a median velocity 
( 1 ) 50 ) more negative than —50 km s -1 . Objects with OH purely 
in emission are thus excluded by this definition, even though one 
of them (PG 1613+658 16 ) shows broad OH profiles extending 
in excess of ~500 km s , suggestive of non-gravitational 


16 FI 345 l-t-1232 also shows broad OH profiles but this may partly be due to 
the binary nature of this object. 


motion . 17 As described in Rupke et al., the —50 km s _1 cutoff 
is used to avoid contamination due to systematic errors and 
measurement errors in wavelength calibration, line fitting (see 
Section 3.3), and redshift determination. The few redshifted 
absorption components detected in our sample have V 50 < 
50-100 km s" 1 (Figure 6 ), suggesting that this cutoff is 
reasonable. 

Using this conservative definition, we detect a molecular wind 
in 26 (70%) of the 37 objects in our sample with measurable 


17 While winds produce blueshifted OH absorption profiles, they are not 
expected to produce shifts of the line profiles if OH is purely in emission, 
unless dust obscuration causes trir 1 . 
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OH 1 19 [i m feature. This fraction is ~65% and 62% if the wind 
velocity threshold is changed to 1)50 = —75 and —100 km s _1 , 
respectively. The requirements to detect the wind in absorption 
against the continuum source and for the absorption velocities 
to be ^—50 km s -1 necessarily bias the sample against winds 
seen edge-on. This wind detection rate is therefore a lower 
limit and our results may be consistent with all ULIRGs having 
such molecular winds. Rupke et al. (2005b) came to a similar 
conclusion regarding neutral winds traced by Na 1 , although the 
sample of Rupke et al. (2005b) is different from the present 
one in that it is made of starburst dominated (U)LIRGs and 
therefore primarily traces starburst-driven winds. The opening 
angle of the molecular outflows inferred for the wind detection 
rate of 70% is ~145°, assuming all objects in our sample have 
an outflow. 

Figure 7 shows no clear dependence of the wind detection rate 
on AGN fraction, AGN luminosity, or star formation luminosity, 
to within the uncertainties of the data (recall that objects with 
OH feature in pure emission are excluded from this analysis; 
these objects are discussed in Section 4.7). In contrast, Rupke 
et al. (2005b) detected Nai winds slightly more frequently in 
ULIRGs (~75%) than in LIRGs (~45%), and attributed this 
difference to higher collimation of the outflows in LIRGs (see 
also Chen et al. 2010). The results of Rupke et al. (2005b) 
refer to starburst dominated LIRGs, which are not adequately 


sampled by our set of objects, so we cannot verify if the depen- 
dence of the wind detection rate on SFR seen in the neutral gas 
is also present in the molecular gas. More relevant is the sample 
of Rupke et al. (2005c), which contains several objects in com- 
mon with the present sample. The lower neutral wind detection 
rate found among Seyfert 2 ULIRGs in Rupke et al. (2005c; 
45% ± 1 1%) is not seen in OH, although the comparison is nec- 
essarily limited to systems in the present sample with q!agn 
90%. We return to the multi-phase nature of these winds in 
Section 5. 

Only four objects in our sample (F 1 2243 — 0036 = N GC 4418, 
15250+3609, F17207-0014, and F22491-1808) show clear 
evidence for inflow based on the detection of a redshifted OH 
absorption feature with median velocity i> 5 o(abs) ^ 50 km s -1 . 
This paucity suggests that the infalling gas is either rare or 
subtends a relatively small fraction of 4 7 r steradians (e.g., 
planar or filamentary geometry). The fact that clear inverted 
P Cygni profiles are not seen in our sample of galaxies adds 
support to this argument. While Gonzalez- Alfonso et al. (2012) 
have shown that the OH 119/rm feature in NGC 4418 is best 
fit as a weak inverted P Cygni profile, this conclusion is driven 
largely by the clear detection of such a profile in [Oi] 63 /um. 
This result serves as a cautionary tale that weak blueshifted OH 
emission may be missed in our single-line fits. Nevertheless, 
this does not change the fact that inverted P Cygni profiles 
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Table 2 

Properties of the OH 119 (tm Profiles 


(Name) 

U 50 (abs) 

U 84 (abs) 

w max (abs) 

Flux a bs 

EQW abs 

V50 (emi) 

«84 (emi) 

FluXemi 

EQW em i 

EQW Tota i 


(km s _1 ) 

(km s -1 ) 

(km s -1 ) 

(Jy km s -1 ) 

(km s -1 ) 

(km s -1 ) 

(km s _1 ) 

(Jy km s _1 ) 

(km s -1 ) 

(km s -1 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

F00509+1225 






171 

363 

150.0 

-65 

-65 

F0 1572+0009 

-724:: 

-892:: 

-1100:: 

-32.8:: 

23:: 

195:: 

345:: 

32.7:: 

-24:: 

0:: 

F05024— 1941 

-183 

-508:: 

-850:: 

-83.6 

83 





83 

F05 189-2524 

-327 

-574 

-850 

-510.0 

59 

189 

351 

424.1 

-51 

10 

07251-0248 

-63 

-255 

-550 

-453.8 

56 





56 

F075 99+65 08 

-459 

-652 

-1000 

-47.7 

35 

315 

808 

74.1 

-53 

-19 

F08572+3915 

-489 

-832 

-1100 

-347.5 

124 





124 

09022-3615 

-153 

-297 

-650 

-265.5 

33 

195 

339 

405.7 

-51 

-17 

F09320+6134 

-9 

-225 

-1200 

-5.0 

83 





83 

F10565+2448 

-267 

-489 

-950 

-2140.7 

156 

99 

273 

1099.2 

-82 

75 

FI 11 19+3257 

-423 

-814: 

-1200: 

-141.6 

124 





124 

Z1 1598—01 14 

-153 

-477 

-800 

-295.2 

121 





121 

FI 2072— 0444 

-69 

-321 

-1200 

-221.9 

135 

243 

393 

137.6 

-87 

51 

F121 12+0305 

-117 

-237 

-400 

-0.8 

52 

243 

447 

0.8 

-48 

2 

F12243— 0036 
F12265+0219 

111 

-21 


-1406.6 

63 





63 

F12540+5708 

-237 

-610 

-1500 

-5510.7 

243 

147 

357 

2956.9 

-134 

113 

13120-5453 

-195 

-520 

-1200 

-5037.9 

113 





113 

F13305— 1739 
F13428+5608 

-201 

-495 

-750 

-679.2 

61 





61 

F13451+1232 






-51 

190 

203.2 

-136 

-136 

F14348— 1447 

-291: 

-508 

-900 

-696.1 

112 

93 

273 

542.1 

-88 

25: 

F14378— 3651 

-219 

-556 

-1200 

-1009.4 

206 

261 

423 

423.6 

-91 

119 

F14394+5332 

F15206+3342 

-291 

-495 

-750 

-106.9 

62 

387 

520 

43.0 

-26 

37 

F15250+3608 

189 

-21 


-491.6 

125 





125 

F15327+2340 

21 

-153 

-700 

-22450.3 

211 





-211 

F15462— 0450 

-225: 

-459: 

-600: 

-161.9 

80: 





80: 

F16504+0228 

-207 

-544 

-1200 

-2126.8 

84 





84 

F17207— 0014 

51 

-165 


-4619.4 
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148 

F19297— 0406 

-231 

-532 

-1000 

-750.5 

119 





119 

19542+1110 

-93: 

-489: 

-700 

-331.9 

69 

243 

441 

191.2 

-41 

29: 

F2055 1-4250 

-381 

-748 

-1200 

-2.1 

70 





70 

F2249 1 — 1808 
F23128— 5919 

99 

3:: 


-0.3 

25: 





25: 

F23233+2817 

-267 

-423 

-500 

-41.3 

30 

201 

339 

79.2 

-58 

-27 

F23365+3604 

-243 

-604 

-1300 

-1348.8 

197 

153 

321 

619.1 

-94 

107 

F23389+0300 
PG 1126-041 

-171 

-285 

-600 

-43.4 

48 

351 

634 

86.5 

-91 

-46 

PG 1351+640 






33 

171 

12.8 

-40 

-40 

PG 1440+356 






-105 

69 

31.1 

-43 

-43 

PG 1613+658 
PG 2130+099 






135 

333 

55.8 

-80 

-80 

Average 

-194 

-444 

-927 

-1605 

103 

202 

401 

439 

-72 

67 

Median 

-204 

-492 

-925 

-340 

84 

198 

354 

197 

-70 

73 


Notes. The uncertainties on the equivalent widths and fluxes (velocities, except u max ) are typically 20% (50 km s -1 ), unless the value is followed by a colon (when the 
uncertainties are 20%-50% and 50-150 km s 1 , respectively) or a double colon (when the uncertainties are larger than 50% and 150 km s -1 , respectively). For u max , 
the typical uncertainty is ±200 km s _1 . Column 1 : galaxy name. Coordinate-based names beginning with F are sources in the IRAS Faint Source Catalog. Column 2: 
t> 5 o(abs) is the median velocity of the fitted absorption profile, i.e., 50% of the absorption takes place at velocities above — more positive than — 1 ) 50 - Column 3: Lotahs) 
is the velocity above which 84% of the absorption takes place. Column 4: i> max (abs) is the maximum extent of the blueshifted wing of the OH 1 19.233 (tm absorption 
profile. Column 5: the total integrated flux for the absorption component(s). Column 6 : the total equivalent width for the absorption component(s). Column 7: U 50 
(emi) is the median velocity of the fitted emission profile. Column 8 : U 84 (emi) is the velocity below which 84% of the emission takes place. Column 9: the total 
integrated flux for the emission component(s). Column 10: the total equivalent width for the emission component(s). Column 11: the total equivalent width for the 
sum of the two components for one line of the OH doublet. 


with strong blueshifted emission are absent from our sample of 
objects (although there are a few known exceptions among other 
ULIRGs, e.g., J. Gracia Carpio et al. 2013, in preparation). This 
means that most of the infalling gas in ULIRGs lies in front 
of the FIR continuum source and therefore does not subtend a 
large solid angle as seen from the source. 


4.4. Kinematics versus Host Galaxy 
Masses and Merger Phases 

The velocities of local neutral gas outflows traced by Nai 
show a tendency to become more negative with increasing 
host galaxy stellar velocity dispersions and stellar masses (e.g.. 
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Figure 3. Histograms showing the overall detection rates of OH 1 19 fim as 
a function of (a) the starburst luminosities, (b) the AGN fraction, and (c) the 
AGN luminosities. The uncertainties on the detection rates assume a binomial 
distribution. No significant trend is found with any of these parameters (the 
statistics are poor at low starburst luminosities). 

Rupke et al. 2002, 2005b, 2005c; Martin 2005; Chen et al. 2010). 
A similar trend is present when using low-ionization tracers in 
Z ~ 1 star-forming galaxies (Erb et al. 2012; cf. Steidel et al. 


2010; Law et al. 2012 at higher redshifts). Figure 8 compares 
all of the velocities of the OH absorption features derived from 
our data (regardless of whether they indicate outflow or inflow) 
with these same host properties. We detect no obvious trend 
with near-infrared derived stellar velocity dispersions or stellar 
masses. Perhaps this is not too surprising since our sample 
covers a considerably narrower range of host galaxy properties 
than the Nai samples. Indeed, the low-z trends between Nai 
outflow velocities and host galaxy properties largely disappear 
when only ULIRGs are considered. Moreover, Rothberg & 
Fischer (2010) and Rothberg et al. (2013) have recently argued 
that the near-infrared derived velocity dispersions of luminous 
IR merger remnants may more closely reflect the properties 
of the young stellar disk than the global properties of the host 
galaxy system as a whole. The noticeably smaller stellar velocity 
dispersions among systems with OH P-Cygni or pure-emission 
profiles is unexpected, and is not seen in the right panels of 
Figure 8 where the host stellar masses are considered. This 
result may be related to the young stellar population bias in 
near-infrared derived velocity dispersions noted by Rothberg & 
Fischer. 

Since nearly all galaxies in our sample are undergoing a 
galaxy interaction, it is natural to ask whether there is a 
dependence of the molecular velocities on merger phase, or 
equivalently, “interaction class” as defined in Veilleux et al. 
(2002). Given the relatively small number of objects with 
detected OH absorption features, we simplified the classification 
into “binaries,” where two distinct nuclei have been detected 
via arcsecond resolution imaging, and “singles,” where the two 
nuclei have coalesced. The results are shown in Figure 9. The 
OH absorption velocities in binaries are not statistically different 
from those in singles. This conclusion is the same whether we 
use i> 5 o(abs), i> 84 (abs), and u max (abs). Recall, however, that IR- 
faint QSOs are not included in this analysis since their OH 
feature is in emission. One of the four IR-faint ULIRGs/QSOs 
with detected OH (PG 1613+658) shows broad OH emission 
profiles and may therefore qualify as having an outflow. We 
return to this point below (Section 4.7). 




Figure 4. Total (absorption + emission) equivalent widths of OH 119 fim (positive values indicate absorption, negative values indicate emission) as a function of (a) 
the AGN fractions and (b) the AGN luminosities. The meanings of the symbols are as follows: filled black squares, filled red triangles, and open red triangles represent 
ULIRGs with OH 119/rm seen purely in absorption, composite absorption/emission, and purely in emission, respectively. The OT2 QSOs, where OH \\9 /im is 
purely in emission, are shown separately as open black circles. The vertical dotted lines indicate the positions of the six sources with undetected OH 119 /im (two of 
these sources have nearly exactly the same AGN fraction). The horizontal gray line marks the null OH equivalent width. 

(A color version of this figure is available in the online journal.) 
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(b) the depths, on a logarithmic scale, of the 9.7 /im features relative to the local continuum measured by Stierwalt et al. (2013; the absolute value of Sgj is actually 
plotted here). In both panels, the strength of the silicate absorption feature increases to the right (see Section 4.2 for definitions of the silicate-related quantities). The 
sign convention for the equivalent widths and the meanings of the symbols are the same as in Figure 4. The vertical dotted lines indicate the positions of the two 
sources with silicate measurements but undetected OH 119 /im. The horizontal gray line marks the null OH equivalent width. 

(A color version of this figure is available in the online journal.) 
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Figure 6. Histograms showing (a) the distributions of the median 50% velocities derived from the multi-Gaussian fits to the OH profiles: absorption components (blue) 
and emission components (red), (b) same as (a) for the 84% velocities, (c) the distributions of the median 50% velocities derived from the absorption components 
in the fits: pure absorption profiles (hatched gray) and P Cygni profiles (dashed blue), (d) same as (c) for the 84% velocities, and (e) the distribution of the terminal 
absorption velocities. Typical uncertainties on U 50 and '-K 4 are ±50 km s -1 and ±200 km s 1 on u max (see Table 2). 

(A color version of this figure is available in the online journal.) 


4.5. Kinematics versus Starburst Luminosities 

Plots of the velocities of local neutral gas outflows traced by 
Nai (on the vertical axis) versus SFRs (on the horizontal axis) 
show an “envelope” with the most negative measured outflow 
velocities increasing with SFRs (e.g., Rupke et al. 2005b; Martin 


2005; Chen et al. 2010). This trend seems driven largely by the 
modest outflow velocities measured in dwarf galaxies (Schwartz 
& Martin 2004) and the large scatter in the outflow velocities 
measured in the ULIRGs (Rupke et al. 2005b; Martin 2005). A 
similar trend may be present at z > 1 when considering wind 
velocities derived from low-ionization species (e.g., Weiner 
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Table 3 

Results from the Statistical Analysis for Sources with t) 5 o(abs) < 50 km s -1 


Parameter 

(i) 

Number of Objects 
(2) 

Break 

(3) 

Pk-s 

(4) 

Ps 

(5) 

(6) 

r 

(7) 

Pr 

(8) 

r 

(9) 

Pr 

(10) 

logf-sB-t^O 

28 

12 

0.11 

0.33 

0.087 

0.23 

0.086 

0.18 

0.37 

logisB-t'84 

28 

12 

0.54 

0.19 

0.33 

0.14 

0.30 

0.19 

0.32 

log LsB-fmax 

28 

12 

0.54 

0.21 

0.27 

0.15 

0.25 

0.23 

0.25 

“AGN-fSO 

28 

50 

0.16 

-0.47 

0.012 

-0.34 

0.012 

-0.47 

0.011 

« AGN-U84 

28 

50 

0.16 

-0.39 

0.038 

-0.26 

0.049 

-0.45 

0.015 

“AGN-Wmax 

28 

50 

0.47 

-0.32 

0.10 

-0.22 

0.097 

-0.35 

0.069 

log Lagn-uso 

28 

11.8 

0.032 

-0.46 

0.013 

-0.34 

0.012 

-0.53 

0.0036 

log LaGN-1’84 

28 

11.8 

0.071 

-0.39 

0.038 

-0.27 

0.041 

-0.50 

0.0070 

log Z-AGN-t'max 

28 

11.8 

0.44 

-0.24 

0.23 

-0.19 

0.16 

-0.33 

0.087 


Notes. Column 1: quantities considered for the statistical test. Column 2: number of objects with no significant inflows (i) 5 o(abs) < 50 km s -1 ). Column 3: value of 
the break used to define the two distributions for the Kolmogorov-Smimov (K-S) test. Column 4: null probability of the K-S test. Column 5: Spearman rank order 
correlation coefficient. Column 6: null probability of the Spearman rank order correlation. Column 7: Kendall’s correlation coefficient. Column 8: null probability of 
Kendall’s correlation. Column 9: Pearson’s linear correlation coefficient. Column 10: Two-tail area probability of Pearson’s linear correlation. Null probabilities less 
than ~0.05 (shown in bold-faced characters) indicate statistically significant trends; those with null probabilities between 0.05 and 0.10 are tentative. 
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Figure 7. Histograms showing the OH wind detection rates, based on the 
presence of an OH 119 /xm absorption feature with a median velocity 1150 more 
blueshifted than —50 km s -1 , as a function of (a) the starburst luminosities, 
(b) the AGN fraction, and (c) the AGN luminosities. The uncertainties on the 
detection rates assume a binomial distribution. No significant trend is found 
with any of these parameters (the statistics are poor at low starburst and AGN 
luminosities). 


et al. 2009), although the situation likely is more complicated 
at these higher redshifts (Steidel et al. 2010; Kornei et al. 2012; 
Martin et al. 2012). Figure 10 shows a similar plot substituting 
the OH velocities for those based on Na I, covering ~2 dex in 
SFRs. Shown on the horizontal axis are the starburst luminosities 
(oc SFR) derived using Psb = (1 — giagn^bol (see Section 2.2). 
No correlation, upper envelope, or increase in scatter within the 


more luminous starbursts is seen in this figure, regardless of 
which velocity (t> 5 o(abs), i> 84 (abs), or v max (abs)) is plotted. 18 If 
anything, there may even be a tendency for the most luminous 
starbursts to have more positive velocities (slower outflows), 
but this trend is not found to be statistically significant: the 
probability, P[null], that the two populations are drawn from 
the same parent population or are uncorrelated is larger than 
~10%, regardless of the statistical test used for the analysis 
(see values listed in Table 3). Recall from Section 2.2 that the 
ULIRGs with the most luminous starbursts in our sample also 
tend to have smaller AGN fractions (<40%). As discussed in the 
next section, these smaller AGN fractions might help explain the 
slower outflows in these objects. Finally, we warn the readers 
that this analysis is necessarily limited by the range in SFRs 
covered by our sample (~2 dex), which is narrower than that in 
low-z Nai D studies. 

4.6. Kinematics versus AGN Fractions 

Figure 11 shows the OH velocities plotted against the 
AGN fractions chagn derived from the /15//30 ratios. A visual 
inspection of this figure suggests that objects with dominant 
AGNs (q!agn ^ 50%) have more negative velocities (faster 
outflows) than objects with dominant starbursts (q!agn ^ 50%). 
However, a K-S test between the velocity distributions of 
dominant AGNs and dominant starburst systems indicates that 
this difference is not statistically significant . 19 The probability 
P[null] that the two populations are drawn from the same parent 
population ~16% for both 1)50 and 1^4 and ~47% for i> max 
(Tables 3 and 4). A similar negative result is found when ccagn 
is varied from 40% to 60% (beyond this range, fluctuations due 
to the small number of objects in one of the two distributions 
dominate the statistics making the results unreliable). However, 
a search for a correlation between the OH kinematics and the 
AGN fraction produces more positive results with P[null] ~ 
l%-5% (not considering the more uncertain u max ; see Table 3). 
The fact that the Pearson’s test for a linear correlation produces 
small P[null] indicates that the data favor a smooth transition 

18 Objects with significant inflows (t^ofabs) ^ 50 km s -1 ) are excluded from 
this analysis, although this conclusion remains the same even when they are 
included. 

15 Once again, objects with significant inflows (loofahs) Js 50 km s -1 ) are 
excluded from this analysis, although this conclusion remains the same even 
when they are included. 
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Figure 8. 50%, 84%, and terminal OH outflow velocities as a function of (a)-(c)-(e) the near-infrared derived stellar velocity dispersions measured by Dasyra et al. 
(2006a, 2006b, 2007; the typical uncertainties on the velocity dispersions are ~50 km s -1 ) and (b)-(d)-(f) the H - band absolute magnitudes of the host galaxies from 
Veilleux et al. (2006) or from Veilleux et al. (2002) assuming H—K' = 0.35 mag (the typical uncertainties on these magnitudes are 0. 1 and 0.3 mag, respectively). The 
meanings of the symbols are the same as in Figure 4. In the lower panels, data points with red circles around them indicate values based on Veilleux et al. (2002). The 
two data points joined by a segment in the left (right) panels correspond to F13451+1232 E and W (F14394+5332 E and W). Typical uncertainties on v$o and i.'s 4 are 
±50 km s _I and ±200 km s 1 on v max . The smaller symbols have larger uncertainties (values followed by double colons in Table 2). 

(A color version of this figure is available in the online journal.) 


Table 4 

Parameters from Linear Regression Analysis for Sources with 
i> 5 o(abs) < 50 km s -1 


X-Y 


Y =aX + b 

a 

b 

« AGN-U50 

- 3.0 ± 1.0 

-92.0 ± 40.3 

«AGN~^84 

- 3.5 ± 1.2 

-332.7 ±61.5 

a AGN~^max a 

- 4.0 ± 1.9 

-732.0 ± 98.4 

log £aGN-E50 

-212 ±72 

2286 ± 838 

log Z.AGN-W84 

-236 ±73 

2310 ±867 

log ^AGN~ '^max 3 

-237 ± 100 

1891 ± 1185 


Notes. X and Y represent the independent and dependent variables, respectively. 
a and b represent the regression coefficient (slope) and regression constant 
(intercept), respectively. For the relationship presented in this table we used 
the ordinary least-square regression of the dependent variable, Y , against the 
independent variable X, OLS(F — X). 

a Note from Table 3 that these correlations are tentative so the parameters listed 
here should be treated with caution. 

with AGN fraction from low-velocity to high-velocity outflows 
rather than a sudden threshold effect. 

The preceding analysis used an equal weighting scheme for 
all data points. The correlations with AGN fractions get slightly 
stronger if we remove from the analysis the velocities with 


the largest uncertainties (those indicated by small symbols in 
Figure 1 1 and double colons in Table 2). 

Note that a positive trend between wind velocities and AGN 
fractions does not necessarily imply that the AGN is the driver 
of these outflows. Once a significant fraction of the obscuring 
material has been swept away, either by an AGN- or starburst- 
driven wind, we expect to be able to see down the core of the 
ULIRG more easily. It therefore becomes easier to observe the 
central high-velocity outflowing material and detect the AGN 
via optical-MIR diagnostics. Orientation effects may also be 
at play in some cases: both the AGN and central high-velocity 
wind are more easily detectable in systems where the dusty 
torus/disk are seen face-on. To test whether the AGN plays a 
role in driving these outflows, we need to compare the outflow 
kinematics with the AGN luminosities; this is done in the next 
section. 

4. 7. Kinematics versus AGN Luminosities 

Figure 12 displays the OH outflow velocities versus the 
AGN luminosities, Lagn = wagn^boi. (Section 2.2). This 
figure shows distinctly more negative velocities in systems 
with larger AGN luminosities. This trend is more evident 
with 1)50 and i>g 4 than with the more uncertain i> max . The null 
probability for a correlation between the OH velocities (vso 
or U84) and log L A gn range from ~0.4% to ~4%, depending 
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Figure 9. 50%, 84%, and terminal OH outflow velocities as a function of 
interaction class: binaries or single sources. The meanings of the symbols 
are the same as in Figure 4. The error bars indicate the lcr scatter in each 
distribution. There are only two ULIRGs with OH purely in emission: one 
of them (F1345 1+1232) is a binary while the other has a single nucleus 
(F00509+1225 = I Zw 1). Typical uncertainties on 1150 and U 84 are ±50 km s -1 
and ±200 km s -1 on u max (see Table 2). 

(A color version of this figure is available in the online journal.) 


on the statistical test (Tables 3 and 4). A visual inspection of 
panels (a) and (b) of Figure 12 also suggests a steepening of the 
relation between outflow velocities and log Lagn fit high AGN 
luminosities: systems with v$o < -400 km s 1 all have AGN 
luminosities larger than log(L^^/L 0 ) = 11.8±0.3. The factor 
of two uncertainty on this luminosity break takes into account 
the measurement uncertainty on this break visually estimated 
from the panels of Figure 12 (estimated to be ±0.2) and the 
uncertainties (±0.3) on the values of ckagn and Lagn for our 
objects. To verify the value of this break more quantitatively, 
K-S tests were carried out between the velocity distributions 
of objects with log (Lagn/Lq) > N and those with log 
(Lagn/Lq) < N, where N was varied from 11.7 to 12.0(beyond 
this range, fluctuations due to the small number of objects in 
one of the two distributions dominate the statistics, making the 
results unreliable). The results from these tests confirm that the 
most significant difference between the two distributions occurs 
when N — 1 1.8 (P[null] ~ 3% for V 50 and 7% for 1^4; Tables 3 
and 4). This break is not statistically significant when using u max . 
These results emphasize the need to use more robust values of 
the outflow velocities, such as V 50 and v$ 4 , when carrying out 
these analyses. 

Once again, the preceding analysis used an equal weighting 
scheme for all data points. The conclusions remain the same 




11 12 13 


log (L /L ) 
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Figure 10. 50%, 84%, and terminal OH outflow velocities as a function of the 
starburst luminosities. The meanings of the symbols are the same as in Figure 4. 
The four objects with significant inflows (50% OH velocities above 50 km s -1 ) 
are not shown here. None of these quantities is significantly correlated with the 
starburst luminosities (P[null] > 0.05; see Table 3). F12072— 0444, indicated 
by an open red triangle, was not included in the evaluations of these correlations 
because it is not clear which of the two nuclei is responsible for the OH 
absorption feature (both nuclei are included in the PACS entrance aperture). 
Typical uncertainties on U 50 and U 84 are ±50 km s -1 and ±200 km s -1 on 
Umax - The smaller symbols have larger uncertainties (values followed by double 
colons in Table 2). 

(A color version of this figure is available in the online journal.) 


if we remove from the analysis the velocities with the largest 
uncertainties (those indicated by small symbols in Figure 12 and 
double colons in Table 2), although in that case the significance 
of the correlations with AGN luminosities become virtually the 
same as with AGN fractions. 

Recall that none of the four ULIRGs/QSOs with AGN 
fraction ctagn 90% and detected OH 119/rm features, 
including all of the OH-detected OT2 IR-faint quasars, shows 
OH absorption components. These objects thus do not appear 
in Figures 10-12 and are excluded from our statistical analysis 
(Table 3). The generally modest widths of the OH emission 
profiles in these objects (PG 161 3+658 is the only exception) 
do, however, suggest that the winds, if present, have modest 
velocities. This result seems inconsistent with the trend of 
increasing outflow velocities with increasing AGN luminosities 
since all four objects with AGN fraction qjagn ~ 90%' have 
AGN luminosities near or above 10 118 L, 0 . Perhaps we are 
seeing the end of the fast wind phase predicted by some models 
(Narayanan et al. 2008; Hopkins et al. 2009). 
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Figure 11. 50%, 84%, and terminal OH outflow velocities as a function of the 
AGN fractions. The meanings of the symbols are the same as in Figure 10. 
The four objects with significant inflows (50% OH velocities above 50 km s -1 ) 
are not shown here. Dotted lines indicate the best linear fits through the data. 
A significant linear correlation is present with V 50 and vs 4 (P[null] 0.05), 
but only tentatively with u max (pfnull] ~0.07; Table 3). The coefficients of 
the correlations are listed in Table 4. FI 2072— 0444, indicated by an open red 
triangle, was not included in the evaluations of these correlations because it is 
not clear which of the two nuclei is responsible for the OH absorption feature 
(both nuclei are included in the PACS entrance aperture). Typical uncertainties 
on V 50 and v %4 are ±50 km s -1 and ±200 km s -1 on u max . The smaller symbols 
have larger uncertainties (values followed by double colons in Table 2). 

(A color version of this figure is available in the online journal.) 



Figure 12. 50%, 84%, and terminal OH outflow velocities as a function of the 
AGN luminosities. The meanings of the symbols are the same as in Figure 10. 
The four objects with significant inflows (50% OH velocities above 50 km s -1 ) 
are not shown here. Dotted lines indicate the best linear fits through the data. 
A significant linear correlation is present with 1)50 and i >84 (/’[null] ^ 0.05), 
but only tentatively with u max (P[null] ~0.09; Table 3). The coefficients of 
the correlations are listed in Table 4. FI 2072— 0444, indicated by an open red 
triangle, was not included in the evaluations of these correlations because it is 
not clear which of the two nuclei is responsible for the OH absorption feature 
(both nuclei are included in the PACS entrance aperture). Typical uncertainties 
on U 50 and v %4 are ±50 km s -1 and ±200 km s -1 on i/ max - The smaller symbols 
have larger uncertainties (values followed by double colons in Table 2). 

(A color version of this figure is available in the online journal.) 


5. DISCUSSION 

Our kinetic analysis of the OH 1 19 fim velocities in nearby 
ULIRGs and QSOs described in Section 4 shows trends of 
increasing molecular wind velocities (becoming more negative) 
with increasing AGN fractions and luminosities. While the 
former may be explained as an obscuration effect where both 
the AGN and central high-velocity outflowing material are more 
easily detectable when the dusty material has been swept away 
or is seen more nearly face-on, the trend with AGN luminosity 
suggests that the AGN is actually playing a role in driving the fast 
winds. 211 Our data favor the existence of a nonlinear (steepening) 
relation between outflow velocities and the logarithms of the 
AGN luminosities. Above 

log(LS/L G ) = 11.8 ±0.3, (3) 


20 We note, however, that for non-spherical geometries, even apparent AGN 
luminosities may be affected by optical depth and projection effects. For 
example, in the simple case of an optically thick ( trir 1 ) torus surrounding 

a MIR-bright compact object, the luminosity can be significantly higher in the 
polar direction. 


the AGN seems to play a dominant role in driving the outflowing 
molecular gas in the objects of our sample. This luminosity break 
coincides approximately with the minimum AGN luminosity 
traditionally used to define a quasar (Mb ^ —23). 

This conclusion only applies to objects with absorption 
signatures of molecular outflows. Fully coalesced ULIRG/QSO 
mergers with AGN fractions q?agn > 90%, including all of the 
OH-detected OT2 IR-faint quasars where OH is seen purely in 
emission, generally show modest line widths and thus modest 
outflows despite their large AGN luminosities. 

The results from S 1 1 on 6 ULIRGs in common with our 
sample indicate that the OH outflows are not only fast but also 
powerful, capable of displacing a significant fraction of the 
entire molecular interstellar medium (ISM) of the host galaxies. 
Given our methods of selection, there is no reason to believe 
that the objects studied here are any different from those in SI 1. 
Since the molecular ISM is the fuel for future star formation 
and galaxy build-up, the quasar-driven outflows detected here 
may therefore have a significant effect on the evolution of their 
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hosts, potentially quenching star formation in their centers. The 
modest OH emission line widths in the IR-faint quasars perhaps 
indicate that AGN feedback subsides once the quasar has poked 
a substantial hole through its natal “cocoon,” although this is 
based on only four objects. 

There has been considerable theoretical support in recent 
years for the idea of quenching star formation in galaxies 
via AGN-driven winds (e.g., Narayanan et al. 2008; Hopkins 
et 2009, and references therein), but relatively little direct 
observational evidence for it, until recently (e.g., Feruglio 
et al. 2010; SI 1; RV11; RV13). This last paper also points 
out the possibility of increasing outflow velocities at higher 
AGN luminosities. Using detailed IFU data on a sample of 
six ULIRGs, all from the QUEST sample, RV13 find that 
fast AGN-driven neutral (Nai D) and/or ionized (Ha) winds 
are only present in systems with log (L^q*/L 0 ) ^ 11.7. 
This AGN luminosity, while only tentatively identified in that 
paper since it is based on a very small sample, is consistent to 
within the uncertainties to that inferred from the larger set of 
Herschel data. This suggests that the conditions to drive fast 
neutral/ionized outflows are approximately the same as for the 
molecular outflows. 

In this context, it is instructive to compare the OH velocities 
measured here with the Na i velocities measured by Rupke et al. 
(2005a, 2005b, and 2005c) and H. B. Krug et al. (2013, in 
preparation) from integrated spectra and RV 1 3 from IFU data 
(Figure 13). Also shown in this figure for comparison are the 
MIR fine-structure line outflow velocities deduced by Spoon & 
Holt (2009) from Spitzer integrated spectra, the Ha velocities 
measured by RV13 from the IFU data, and the Hi 21 cm 
velocities measured by Teng et al. (2013) from GBT integrated 
spectra. While the various velocities are measured over different 
scale (<200 pc for the OH velocities according to Sll and 
<1—2 kpc for the neutral and ionized gas velocities), general 
trends are observed: the outflow velocities of the molecular gas 
are often comparable to the velocities of the neutral gas, but 
they are generally smaller (more positive) than the velocities of 
the ionized gas (panels (b) and (c) in Figure 13). The similarity 
between the OH and Nai velocities is particularly evident in 
panel (b) of Figure 1 3 where the sometimes very large nuclear 
Nai and Ha outflow velocities, produced on sub-parsec scales 
(e.g., Mrk 23 1 ; Veilleux et al. 2013; Figure 13(a)), are excluded 
from the means. This panel also reiterates the result found by 
RV13 that the ionized gas is often, but not always, the fastest 
component of the outflow. While dust obscuration may explain 
some of these variations (especially for the highly ionized gas 
found near the center of AGN-dominated ULIRGs; e.g., blue 
wings in [Nev] emission profiles of Spoon & Holt 2009), it 
likely does not explain all of the spatial variations found in the 
IFU data. 

These multi-phase comparisons may help constrain the pro- 
cesses involved in accelerating the material entrained in the 
wind. In most cases, the acceleration a — F /m is expected to 
depend strongly on the ISM phase. For instance, in the case 
of ram pressure acceleration of clouds in a fast wind or a jet 
(e.g., StreTnitskii & Sunyaev 1973; Chevalier & Clegg 1985), 
the force on a spherical cloud of mass M c i oud and density p c i oud 
is expected to scale as the area of the cloud or (M c i ou d/p c ioud) 2 ^ 3 
and thus we expect the acceleration to scale as (M c i ou dp 2 loud ) _1 ' /3 
(we neglect gravity). Naively, one would expect the molecular 
clouds to be considerably denser and more massive, on average, 
than the ionized gas clouds so one would expect the molecu- 
lar material to always move at lower velocities. This simple- 


minded argument may be relaxed if each accelerated cloud has 
a multi-phase structure (e.g., ionized skin or halo surrounding a 
neutral/molecular core) and thus contributes all at once to the 
molecular, neutral, and ionized phases of the outflow. But the 
significant difference in velocities often seen between the ion- 
ized and molecular gas phases (Figure 13) indicates that this sit- 
uation cannot apply to the bulk of the outflowing material since 
the difference in velocities between the various phases would 
imply severe shearing/shredding and very short cloud lifetimes 
(~Rcioud/Au ~ 10 4 (R c i ou d/pc)(100 km s _1 /Ai0 years, where 
Rcioud is the cloud radius and An is the difference in velocities 
between the various gas phases). 

Another possible source of acceleration for these winds is 
UV-IR radiation pressure. It has the distinct characteristic of 
being less efficient at accelerating ionized material than dusty 
neutral/molecular material due to the different sources of 
opacity. Within the dust sublimation radius, R su bi ~ L 4 g pc 
(where a dust sublimation temperature T su bi — 1200 K was 
assumed), the dominant source of opacity is electron scattering 
and the corresponding luminosity condition to drive a (spherical) 
wind is the well-known Eddington luminosity. Estimates of the 
BH masses Mbh in our sample galaxies range from ~0.5-5 x 
10 8 M q (e.g., Vestergaard & Peterson 2006; Veilleux et al. 
2009a), corresponding to Eddington luminosities Lem — 1.3 x 
10 46 (Mbh/10 8 ) = 0.65-6.5 x 10 46 erg s -1 . Systems near 
the AGN luminosity break identified in our data are therefore 
accreting at sub-Eddington values T = ~ 0.04-0.4. 

Outside the dust sublimation radius, the much larger dust 
opacity prevails and UV-IR radiation pressure from the AGN 
may drive significant winds. In this situation where dust opacity 
dominates, Murray et al. (2005) have derived a handy expres- 
sion for the limiting, Eddington-like luminosity above which 
momentum deposition from the quasar (and/or starburst) is pre- 
sumed to be enough to clear most of the gas out of galactic nuclei 
(i.e., this is a criterion to “blow away” the gas rather than just a 
“blow out” in the nomenclature of De Young & Heckman 1994): 

L m - (4f g c)/G)a 4 

~ 3 x 10 46 (/ g /0.1)(cr/200 km s -1 ) 4 erg s -1 , (4) 

where / s is the gas fraction of the host and a is the stellar velocity 
dispersion, a measure of the depth of the gravitational potential 
well. In our sample galaxies, f g ~ 0.1 on average, although 
there is considerable scatter among ULIRGs (e.g., Downes & 
Solomon 1998; Combes et al. 2013, and references therein), 
while a ranges from ~ 120 to ~ 280 km s~ 1 with an average value 
of ~200 km s -1 (see Figure 8; as mentioned in Section 4.4, this 
value can be considered a lower limit; Rothberg & Fischer 2010; 
Rothberg et al. 2013). The steep dependence of Lm on a implies 
that L^gn — 2%-100% of L M , with an average value of ~10%, 
i.e., the AGN luminosity break above which the AGN plays a 
dominant role in driving the outflowing molecular gas in local 
ULIRGs is generally lower than the Eddington-like luminosity 
derived by Murray et al. (2005). 

The fact that L^gn is lower than Lm is not surprising since 
the latter quantity is the maximum luminosity an object may 
have before blowing away all of the gas from a galaxy via ra- 
diation pressure on dust. In addition, the theoretically derived 
Lm is subject to a number of issues that may overestimate its 
value in a real system: (1) First, it is important to note that 
the expression for Lm is an order-of-magnitude estimate of the 
actual value. A number of simplifying assumptions were made 
to derive this expression, e.g., the ISM is approximated as an 
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Figure 13. Multi-phase comparison of the outflow velocities. The 50% OH outflow velocities are compared with (a) the 50% Nai D outflow velocities measured by 
Rupke et al. (2005a, 2005b, 2005c) and H. B. Krug et al. (2013, in preparation) from integrated spectra (in F05024— 1941, F07599+6508, and 12540+5708, there are 
multiple Nai components seen at different velocities), (b) the spatially averaged 50% Nai D (encircled) and Ha outflow velocities measured by RV13 outside of the 
nucleus using the IFU data, (c) the outflow velocities measured by Spoon & Holt (2009) from the emission profiles of [Nem] 15.5 fim (encircled) and [Ne v] 14.3 fim 
in integrated Spitzer spectra, and (d) the H i 21 cm outflow velocities measured by Teng et al. (2013) from GBT data. The meanings of the symbols are the same as in 
Figure 4. The red dotted line is the 1 : 1 line. 

(A color version of this figure is available in the online journal.) 


isothermal sphere of gas with constant gas fraction. The ISM in 
ULIRGs is undoubtedly much more complex. Disky ISM distri- 
butions, such as those presumed to exist in a number of ULIRGs 
(e.g., Downes & Solomon 1998), offer less resistance to winds 
along the disk rotation axis for a given Lagn- A good exam- 
ple is the wide-angle galactic wind in Mrk 231 (RV11, RV13), 
which is inferred to be oriented perpendicular to the central 
molecular disk, yet it entrains a substantial fraction of the gas 
in the nuclear region (i.e., it is closer to a “blow away” event 
than just a “blow out,” so Lm is relevant here). In situ forma- 
tion of the accelerated clouds (Faucher-Giguere et al. 2012) 
and/or energy-conserving acceleration due to inefficient 
cooling in the shocked winds of the AGN (Faucher-Giguere 
& Quataert 2012) may also effectively reduce the value of Lm- 
(2) The expression for the maximum luminosity only takes into 
account the radiation pressure from the AGN. Momentum injec- 
tion from loosely collimated jets (e.g., Teng et al. 2013) would 
reduce the AGN luminosity needed for efficient molecular out- 
flow. (3) While not dominant in systems with AGN luminosities 
above the powerful starbursts known to exist in these 

systems also contribute energy and momentum to driving the 
winds (including possibly significant contributions from cos- 
mic rays). We have neglected the starburst contributions in the 
present treatment. This would relax the requirements on the 
AGN, effectively reducing Lm from the AGN by a factor of up 
to ~2. 

Observationally, the gas mass fraction f g is difficult to 
determine in ULIRGs and depends on the uncertain CO-to- 
FU conversion factor. The value of f g likely varies considerably 
from one galaxy to the other in our sample. This is particularly 
true of the IR-faint QSOs, although they are still undoubtedly 


gas-rich (e.g., Scoville et al. 2003). Given the uncertainties 
on f g and cr, the large galaxy-to-galaxy variations in these 
quantities, and the steep dependence of Lm on a, it is in fact 
not surprising to find that our data (Figure 12) favor a trend, 
albeit nonlinear, between OH velocities and AGN luminosities 
rather than a clear-cut AGN luminosity threshold for powerful 
molecular outflows. Moreover, the exact shape of this outflow 
velocity-AGN luminosity relation likely applies only to local 
gas-rich ULIRGs. Comparisons with other types of objects with 
different f g , cr, and redshifts have to be done with great care. The 
gas mass fraction of local giant spirals is ~7%-10% (Leroy et al. 
2008; Saintonge et al. 2011), similar to that of local ULIRGs, 
but there is mounting evidence that it increases with redshifts, 
reaching ~34% ± 5% and 44% ± 6% among “main-sequence” 
galaxies at z ~ 1.2 and ~2.3, respectively (Tacconi et al. 2010; 
Daddi et al. 2010). New observations also indicate that the gas 
mass fraction increases with redshift among ULIRGs, by a factor 
of 3 ± 1 from z = 0 to 1 (Combes et al. 2013). These results 
fit well with the large gas mass fractions (~40%) previously 
derived in z ~ 2-3.4 submillimeter galaxies (Tacconi et al. 
2006). For objects with the same stellar velocity dispersion, L M 
at high redshifts is thus expected to be three to four times larger 
than in local objects. Future ALMA observations should be able 
to verify this assertion. 

6. CONCLUSIONS 

Herschel / PACS data from the SHINING guaranteed-time key 
project were combined with data from cycles 1 and 2 open-time 
programs to carry out a systematic search for molecular outflows 
in a sample of 43 nearby ULIRGs and QSOs using OH 1 19 fim. 
The main results from this surveys are: 
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1. The OH 119 j/m feature was detected in 37 (86%) of the 
43 sample galaxies. 

2. The nature of the OH 119 /tm feature (emission, absorption, 
or both) depends on the properties of the galaxies: OH 
emission is stronger relative to OH absorption in quasar- 
dominated systems, becoming completely dominant (i.e., 
pure OH emission) in objects where the quasar contributes 
more than ~90% to the bolometric luminosity (called the 
AGN fraction for short). This behavior is similar to that of 
the silicate 9.7 p , m feature studied with Spitzer and the H i 
21 cm feature measured with GBT in a subset of ULIRGs 
and QSOs. 

3. Unambiguous evidence for molecular outflows, based on 
the presence of absorption features with median velocities 
more blueshifted than —50 km s -1 , is seen in 26 (70%) of 
the 37 targets with detected OH 119/rm. Given possible 
projection effects, this wind detection rate is a lower limit 
and is consistent with all ULIRGs having molecular winds 
with an average opening angle ~ 145°. In contrast, absorp- 
tion features with median velocities above +50 km s -1 , 
indicative of inflow, are detected in only four objects, none 
of which show clear inverted P Cygni profiles. This result 
suggests that molecular inflows are either rare or subtend 
a relatively small fraction of 4 it steradians (e.g., planar or 
filamentary geometry). 

4. Typical median outflow velocities are 200 km s -1 , but 
terminal outflow velocities of - — 1000 km s -1 are detected 
in several objects. Deprojected outflow velocities will be 
even larger. The kinematics of these molecular outflows 
do not show any obvious dependence on the properties of 
the host galaxies (e.g., near-infrared derived stellar velocity 
dispersions, stellar masses) and SFRs, although our sample 
covers a range of properties that is narrower than that of 
low-z non-molecular wind studies, where trends have been 
observed with these quantities. 

5. There is a tendency for objects with dominant and luminous 
AGN to show more blueshifted OH velocities (faster 
outflows). Our data favor a steepening of the relation 
between outflow velocities and the logarithms of the AGN 
luminosities above log (Lagn/^o) = H-8 ± 0.3, although 
better statistics are needed, particularly at lower starburst 
and AGN luminosities, to confirm this nonlinearity. While 
the trend between outflow velocities and AGN fractions 
may be explained as an obscuration effect, where both 
the AGN and central high-velocity outflowing material are 
more easily detectable when the dusty material has been 
swept away or is seen more nearly face-on, the dependence 
of the velocities on the AGN luminosities suggests that the 
AGN is playing a role in driving the fast winds. Our results 
emphasize that the use of the terminal velocities as the sole 
measures of outflow velocities should be avoided due to the 
much larger uncertainties on these measurements. 

6. The OH emission profiles in three of the four late-stage 
ULIRG/QSO mergers with AGN fractions above ~90% are 
narrow, despite their large AGN luminosities, suggesting 
that the winds in these objects, if present, often have 
modest velocities. Quasar feedback thus seems to subside 
among mergers that have cleared a path through the dusty 
circumnuclear material. 

We thank the referee for suggesting changes that helped im- 
prove this paper. Support for this work was provided by NASA 
through Herschel contracts 1427277 and 1454738 (S.V. and 


M.M.) and contracts 1364043, 1435724, and 1456609 (J.F.). 
S.V. also acknowledges support from the Alexander von Hum- 
boldt Foundation for a “renewed visit” to Germany following 
up the original 2009 award, and thanks the host institution, 
MPE Garching, where a portion of this paper was written. 
E.G. -A. is a Research Associate at the Harvard-Smithsonian 
Center for Astrophysics, and thanks the support by the Span- 
ish Ministerio de Economfa y Competitividad under projects 
AYA20 10-21697-C05-0 and FIS2012-39162-C06-01. Basic re- 
search in 1R astronomy at the Naval Research Laboratory is 
funded by the U.S. Office of Naval Research. This research 
made use of PySpecKit, an open-source spectroscopic toolkit 
hosted at http://pyspeckit.bitbucket.org. This work has made 
use of NASA’s Astrophysics Data System Abstract Service and 
the NASA/IPAC Extragalactic Database (NED), which is op- 
erated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. 

REFERENCES 

Aalto, S., Garcia-Burillo, S., Muller, S., et al. 2012, A&A, 537, 44 

Chen, Y.-M., Tremonti, C. A., Heckman, T. M., et al. 2010, AJ, 140, 445 

Chevalier, R. A., & Clegg, A. W. 1985, Natur, 317, 44 

Cicone, C., Feruglio, C., Maiolino, R., et al. 2012, A&A, 543, 99 

Cicone, C., Maiolino, R., Sturm, E., et al. 2013, A&A, submitted 

Combes, F., Garcia-Burillo, S., Braine, J., et al. 2013, A&A, 550, 41 

Daddi, E., Bournaud, F., Walter, F., et al. 2010, ApJ, 713, 686 

Dasyra, K. M., Tacconi, L. J., Davies, R. 1., et al. 2006a, ApJ, 638, 745 

Dasyra, K. M., Tacconi, L. J., Davies, R. I., et al. 2006b, ApJ, 651. 835 

Dasyra, K. M., Tacconi, L. J., Davies, R. I., et al. 2007, ApJ, 657, 102 

De Young, D„ & Heckman, T. M. 1994, ApJ, 431, 598 

Downes, D., & Solomon, R M. 1998, ApJ, 507, 615 

Erb, D. K„ Quider, A. M„ Henry, A. L„ & Martin, C. L. 2012, ApJ, 759, 26 

Fabian, A. C. 1999, MNRAS, 308, L39 

Faucher-Giguere, C.-A., & Quataert, E. 2012, MNRAS, 425, 605 
Faucher-Giguere, C.-A., Quataert, E„ & Murray, N. 2012, MNRAS, 420, 1347 
Feruglio, C., Maiolino, R., Piconcelli, E., et al. 2010, A&A, 518, L155 
Fischer, J., Sturm, E., Gonzalez-Alfonso, E., et al. 2010, A&A, 518, L41 (F 1 0) 
Ginsburg, A., & Mirocha, J. 2011, Astrophysics Source Code Library, record 
ascl: 1109.001 

Gonzalez-Alfonso, E., Fischer, J., Gracia-Carpio, J., et al. 2012, A&A, 541, 4 
Gonzalez-Alfonso, E.. Fischer, J., Gracia-Carpio, J., et al. 2013, A&A, submitted 
Hailey-Dunsheath, S., Sturm, E., Fischer, J., et al. 2012, ApJ, 755, 57 
Heckman, T. 2002, in ASP Conf. Proc. 254, Extragalactic Gas at Low Redshift, 
ed. J. S. Mulchaey & J. Stocke (San Francisco, CA: ASP), 292 
Hopkins, P. F., Murray, N., & Thompson, T. A. 2009, MNRAS, 398, 303 
Iwasawa, K., Sanders, D. B„ Teng, S. H., et al. 2011, A&A, 529, 106 
Kauffmann, G., Heckman, T. M„ White, S. D. M„ et al. 2003, MNRAS, 
341,54 

King, A. 2003, ApJL, 596, L27 

Kornei, K. A., Shapley, A. E., Martin, C. L., et al. 2012, ApJ, 758, 135 
Krug, H. B„ Rupke, D. S. N., & Veilleux, S. 2010, ApJ, 708, 1 145 
Laurent, O., Mirabel, I. F., Charmandaris, V., et al. 2000, A&A, 359, 887 
Law, D. R., Steidel, C. C., Shapley, A. E., et al. 2012, ApJ, 745, 85 
Leroy, A. K„ Walter, F„ Brinks, E„ et al. 2008, AJ, 136, 2782 
Markwardt, C. B. 2009, in ASP Conf. Ser., 411, ADASS XVIII, ed. D. A. 

Bohlender, D. Durand, & P. Dowler (San Francisco, CA: ASP), 251 
Martin, C. L. 2005, ApJ, 621, 227 
Martin, C. L. 2006, ApJ, 647, 222 

Martin, C. L.. Shapley, A. E., Coil, A. L., et al. 2012, ApJ, 760, 127 
Murray, N., Quataert, E., & Thompson, T. A. 2005, ApJ, 618, 569 
Narayanan, D., Cox, T. J., Kelly, B., et al. 2008, ApJS, 176, 331 
Netzer, H., Lutz, D., Schweitzer, M., et al. 2007, ApJ, 666, 806 
Newman, S., Genzel, R„ Forster-Schreiber, N. M., et al. 2012, ApJ, 761, 43 
Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, LI 
Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2 
Rothberg, B., & Fischer, J. 2010, ApJ, 712, 318 

Rothberg, B., Fischer, J., Rodrigues, M., & Sanders, D. B. 2013, ApJ, 767, 72 
Rubin, K. H. R„ Weiner, B. J., Koo, D. C, et al. 2010, ApJ, 719, 1503 
Rupke, D. S. N„ & Veilleux, S. 2011, ApJL, 729, L27 (RV11) 

Rupke, D. S. N„ & Veilleux, S. 2013a, ApJ, 768, 75 (RV13) 

Rupke, D. S. N„ & Veilleux, S. 2013b, ApJL, 775, L15 


20 


The Astrophysical Journal, 776:27 (21pp), 2013 October 10 


Veilleux et al. 


Rupke, D. S., Veilleux, S., & Sanders, D. B. 2002, ApJ, 570, 588 

Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005a, ApJS, 160, 87 

Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005b, ApJS, 160, 115 

Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005c, ApJ, 632, 751 

Saintonge, A., Kauffmann, G., Wang, J., et al. 2011, MNRAS, 415, 61 
Sanders, D. B„ & Mirabel, I. F. 1996, ARA&A, 34, 749 
Sanders, D. B., Soifer, B. T., Elias, J. H., et al. 1988, ApJ, 325, 74 
Schwartz, C., & Martin, C. L. 2004, ApJ, 610, 201 
Schweitzer, M., Groves, B., Netzer, H., et al. 2008, ApJ, 679, 101 
Scoville, N. Z., Frayer, D. T., Schinnerer, E., & Christopher, M. 2003, ApJL, 
585, L105 

Spinoglio, L., Malkan, M. A., Smith, H. A., Gonzalez-Alfonso, E., & Fischer, 
J. 2005, ApJ, 623, 123 
Spoon, H„ & Holt, J. 2009, ApJL, 702, L42 
Steidel, C. C., Erb, D. K„ Shapley, A. E., et al. 2010, ApJ, 717, 289 
Stierwalt, S., Arrnus, L., Surace, J. A., et al. 2013, ApJS, 206, 1 


Strauss, M„ Huchra, J. P., Davis, M., et al. 1992, ApJS, 83, 29 

Strel'nitskii, V. S., & Sunyaev, R. A. 1973, SvA, 16, 579 

Sturm, E., Gonzalez-Alfonso, E., Veilleux, S., et al. 201 1, ApJL, 733, L16 (SI 1) 

Tacconi, L. J., Genzel, R., Neri, R., et al. 2010, Natur, 463, 781 

Tacconi, L. J., Neri, R., Chapman, S. C., et al. 2006, ApJ, 640, 228 

Teng, S. H„ Veilleux, S., & Baker, A. J. 2013, ApJ, 765, 95 

Veilleux, S. 2012, JPhCS, 372, 012001 

Veilleux, S„ Cecil, G., & Bland-Hawthom, J. 2005, ARA&A, 43, 769 
Veilleux, S., Kim, D.-C., Peng, C. Y„ et al. 2006, ApJ, 643, 707 
Veilleux, S„ Kim, D.-C., Rupke, D. S. N„ et al. 2009a, ApJ, 701, 587 
Veilleux, S„ Kim, D.-C., & Sanders, D. B. 1999, ApJ, 522, 113 
Veilleux, S., Kim, D.-C., & Sanders, D. B. 2002, ApJS, 143, 315 
Veilleux, S„ Rupke, D. S. N„ Kim, D.-C., et al. 2009b, ApJS, 182, 628 (V09) 
Veilleux, S., Trippe, M., Hamann, F., et al. 2013, ApJ, 764, 15 
Vestergaard, M., & Peterson, B. M. 2006, ApJ, 641, 689 
Weiner, B„ Coil, A. L.. Prochaska, J. X., et al. 2009, ApJ, 692, 187 


21 


